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Abstract 

The way a charm-quark fragments into a charmed hadron is a challenging problem both for 
the theoretical and the experimental particle physics. Moreover, in neutrino induced charm- 
production, peculiar processes occur such as quasi-elastic and diffractive charm-production 
which make the results from other experiments not directly comparable. We present here a 
method to extract the charmed fractions in neutrino induced events by using results from e + e~, 
ttN, 7 TV experiments while taking into account the peculiarities of charm-production in neutrino 
interactions. As results, we predict the fragmentation functions as a function of the neutrino 
energy and the semi-muonic branching ratio, B^, and compare them with the available data. 



1 Introduction 



The problem of how a charm-quark fragments into a charmed hadron is challenging both from the theoretical 
and the experimental point of view. Indeed, perturbative QCD is not applicable at energies comparable with 
the charm-quark mass. Therefore, parameterizations have to be used and the parameters are determined 
experimentally. In the experiments the direct identification of the charmed hadron in the final state is 
only possible through the visual observation of the hadron decay and the measurement of the kinematical 
variables, practically feasible only with nuclear emulsions. 

In the following, we focus on the so-called charm-production fractions (//j's); i.e. the probability that a 
charm-quark fragments into a charmed hadron h (= D , D + , D s , A c ). 

In this paper we review all existing data on charm-production fractions as measured by e + e~, ttN and 
"fN experiments and predict fh in neutrino induced charm-production. Indeed, data on charmed fractions 
in interactions induced by neutrinos are rather scarce. Only one experiment, E531 M, measured fh with a 
statistics of 122 events with an identified charmed hadron in the final state. 

After a thorough overview of all available data, we present a method to extract fh for neutrino exper- 
iments. The major difference of this method with respect to other presented in the past, i.e. see Ref. ||, 
is that/Vs from e + e — are not extrapolated to neutrino experiments straightway, but some peculiarities of 
z^-induced interactions are accounted for. Indeed, neutrinos may undergo to quasi-elastic and diffractive 
processes with production of charmed hadrons in the final state. These results are also used to estimate 
the semi-muonic branching ration of charmed hadrons. Finally, we discuss the statistical and systematic 
uncertainties associated with our predictions and compare them with the E531 results. 

2 Charm-production fraction in DIS interactions 

2.1 Measurement of the D + / D° ratio 

In this Section we give an overview of the available data on D + /D° measurements in e + e _ , irN, pN and 
jN experiments. Furthermore, we also study the correlated variable Fy, which is defined as V/ (V + P) , 
where V and P signify vector and pseudo-scalar charmed mesons, respectively. When needed, we recomputed 
Fy by using the latest charmed- hadron branching-ratios (BR). Finally, we extract for the first time Fy from 
neutrino induced charm-production data. 

In the following we assume that D*° and D* + production cross-sections are equal at the parton level, 
as well as the direct production of D° and D + . 

From here on, we indicate with D the sum of both prompt and decayed (D* — > D) D meson production. 

Under the previous assumptions, the measurement of D + / D° and D* j D ratios allows us to extract 
Fy. By using the formulae given in Ref. 0, Fy can be extracted from the relations 
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6.5 ±0.8 
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9.3 ±2.4 


1.9 ±0.9 










4.16 [: 




7.8 ±0.8 


2.1 ±0.7 










5.20 [ 




4.7 ±0.8 


1.7 ±0.4 










10.55 




1.36 ±0.16 


0.57 ±0.08 


0.78 ±0.17 


0.65 ±0.05 


0.41 ±0.13 


0.20 ±0.08 


10.55 




1.14±0.15 


0.56 ±0.08 




0.54 ±0.08 


0.42 ±0.07 


0.45 ±0.07 



Table 1: Measured cross-sections for the production of charmed hadrons at center of mass energies 
in the range 4.03 4- 10.55 GeV. 

where R x = D+/D , R 2 = (D° from D*+)/D° and = BR{D* + D°ir+) = (67.7 ± 0.5)% @. 

2.1.1 e + e~ experiments 

The basic principle to reconstruct e + e _ events with charmed hadrons in the final state is common to all the 
experiments. Being a cc pair produced in the annihilation, one charmed hadron is used to tag the events, 
while the other one to study the decay properties. The decay modes used to tag the event are [] 

D *+ _> D o n + _^ (K- 7r +)vr+ (0.026 ± 0.006) 



D° R-TT+ (0.0383 ± 0.0009) 
D° -> K^TT+vr+Tr" (0.0749 ± 0.0031) 



D + -> K"7r + 7r + (0.090 ± 0.006) 
where the corresponding branching ratios are also given in brackets H . 

Results on the total cross-sections for inclusive production of the charmed particles D*°, D* + , D° and 
D + at various \/s are shown in Table [l]. 

A complete review of the probabilities (/(c — > C)) that a c-quark fragments into a D*, D°, D + and 
other charmed hadrons as measured in Z° decays is given in Ref. § and reported in Table H. 

From Tables and || we can extract both R\ and Fy, the latter being estimated by using Eqs. ([!]) 
and (||). The results are given in Table ^ and show that within the experimental errors, both R\ and Fy 
are independent of the energy. 



2.1.2 irN and ^yN experiments 

Several experiments have studied charm-production and extracted R\ and Fy by using 7T beams of different 
energies impinging onto different targets. A non-exhaustive list of all available data in 7T./V experiments is 
given in Table ^. 

1 Here and in the following by D + (D* + ) we implicitly indicate also the charge conjugated decay modes. 
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ALEPH 


DELPHI 


OPAL 


Weighted average 


/(c-B u )(%) 


55.9 ±1.7 ±1.5 


54.5 ±1.5 ±3.2 


58.8 ±4.1 ±4.0 


55.8 ±1.8 


/(c- £>+)(%) 


23.8 ± 0.8 ± 1.3 


22.6 ± 0.8 ± 1.4 


23.1 ±3.0 ±2.0 


23.2 ± 1.0 


/(c -!>,)(%) 


11.5±1.9±0.7 


12.4± 1.1± 1.2 


9.0±2.4±1.1 


11.5 ± 1.1 


/(c-A c )(%) 


7.8 ±0.8 ±0.4 


8.6 ± 1.8 ± 1.0 


4.8 ±2.2 ±0.8 


7.6 ±0.8 


/(c -£>**)(%) 


23.3 ± 1.0 ±0.8 


25.5 ±1.5 ±0.6 


22.8 ±0.9 


23.4 ±0.7 



Table 2: Measured and averaged probabilities that a charm-quark fragments into D°, D + , D s , A c 
and D* + in e + e~ annihilation at -y/s = M z o. 



(GeV) 


Ri 


R 2 


Fv(R%) 


F V {R 2 ) 


4.03 [ 




0.33 ±0.06 




0.74 ±0.10 




4.14 [ 




0.20 ±0.11 




0.91 ±0.23 




4.16 [ 




0.27 ±0.09 




0.85 ±0.16 




5.20 [ 




0.36 ±0.10 




0.70 ±0.16 




10.55 [ 


fs| 


0.45 ±0.06 


0.32 ±0.04 


0.56 ±0.08 


0.70 ±0.13 


M z o 


0.42 ±0.03 


0.28 ±0.01 


0.60 ±0.04 


0.57 ±0.03 


Weighted average 


0.39 ±0.02 


0.29 ±0.01 


0.65 ±0.03 


0.60 ±0.03 



Table 3: Measured Ri, R 2 , Fy in e + e experiments as a function of \fs and their averaged values. 

The NA14/2 photo-production experiment j| measured both R 1 = 0.37 ± 0.10 and R 2 = 0.26 ± 0.04 
from which we can extract the weighted average Fy = 0.57 ± 0.09. 

From these data we can conclude that, within the experimental errors, Ri is both process- and energy- 
independent. 

2.1.3 uN experiments 

Recently two measurements which allowed us to extract for the first time Fy from neutrino experiments 
became available. In Ref. jl2| the CHORUS Collaboration presented a measurement of the production 
rate of D° based on a sample of about 26000 charged-current events interactions located and analyzed 
so far in the target emulsions. After reconstruction of the event topology in the vertex region, 283 D° 
decays were observed with an estimated background of 9.2 events from K° and A decays. The CHORUS 
Collaboration measured the D° production cross-section times BR(D° — ► V2) + BR(D° — ► VA)] p|. 
The total cross-section has been extracted by accounting for the D° decays into all neutrals. The value we 
used is BR(D° -> all neutral) = (25 ± 5)% @. 

Therefore, the D° production cross-section normalized to Vu charged-current (CC) interactions is 

^-1 = (2.65 ±0.18 ±0.24 ±0.50) x 10~ 2 (3) 

at 27 GeV average energy. Notice that this measurement includes both D° prompt and D° from 
the decay of D* mesons. 

The D* + production in ZAj charged-current interactions has been measured, with a similar ZAj beam, 
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Ri 


R2 


Fv(Ri) 


F v (Ri) 


WA92: 350 GeV/c ir- 
on Cu, W 


0.423 ± 0.012 


0.280 ±0.015 


0.60 ± 0.02 


0.57 ± 0.04 


E769: 250 GeV/c ir- 
on Be, Al, Cu, W 


0.419 ±0.043 


0.222 ±0.031 


0.60 ±0.06 


0.42 ± 0.08 


E769: 210 GeV/c ir- 
on Be, Al, Cu, W 


0.258 ±0.058 




0.87 ±0.11 




E653: 600 GeV/c ir- 
on emulsion 


0.393 ±0.032 




0.64 ±0.05 




NA32: 230 GeV/c ir- 
on Cu 


0.422 ±0.033 


0.262 ±0.026 


0.60 ±0.05 


0.52 ±0.07 


NA32: 200 GeV/c ir- 
on Si 


n znq+ - 123 

u.4«jy_ 094 


0.319 ±0.095 


0.58 ±0.16 


0.69 ±0.30 


NA27: 360 GeV/c ir- 
on H 


0.564 ±0.171 




0.41 ±0.21 




Weighted average 


0.415 ±0.010 


0.268 ±0.012 


0.611 ±0.015 


0.541 ± 0.033 



Table 4: Summary of available R±, R2 and Fy measurements extracted from 7rN experiments. 

by the BEBC |J] and NOMAD @j experiments to be (1.22 ± 0.25) x 10~ 2 and (0.79 ± 0.20) x 10~ 2 , 
respectively. The weighted average D* + production rate normalized to charged-current interactions is 

— > - = (0.96 ± 0.16) x 10" 2 (4) 

&CC 

From the measured ratios (||) and (|J), and by knowing B+, we can compute R2 = 0.25 ± 0.06. From 
the latter value and from Eq. (0) we can extract 



F v = 0.50 ± 0.12 

It is worth to notice that R± and Fy extracted from neutrino experiments can be compared straightway 
to e + e~, irN and results, being D + and D° either produced promptly or from the decay of prompt 
D*+ and D*°. Namely, processes peculiar of v interactions do not affect R\ and Fy 0. 



2.1.4 Summary and discussion of all available data on R\ and Fy 

From results reported in the previous Sections, we can argue that within the experimental errors R\ is 
constant over a wide range of energies (y/s ~ 4 ± 90 GeV) and independent of the process. The constant 
behavior of R\ down to yfs ~ 4 GeV can be derived with simple arguments: 

• the masses of D+(1869.3), D°(1864.5), L>* + (2010.0) and L>*°(2006.7) are very similar. Therefore, 
the threshold suppression of D* mesons which tends to enhance D + contribution, is very little; 

• whatever D* / D meson is produced a pion should be always created. Therefore, all charmed mesons 
have the same threshold behavior, which cancels out in the ratio. 

2 In neutrino interactions may also be produced diffractively but, due to the V c d suppression, its rate is 

expected to be about (1.6 ± 0.3) x 10 -4 with respect to CC interactions and therefore negligible. 
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Fy (mcas) 


Fy (UCLA) 


Fy(JETSET) 


Fy (HERWIG) 




= 4.03 GcV 


0.74 ±0.10 


0.47 


0.52 


0.54 


e + e~ 


V^ = 4.14 GeV 


0.91 ±0.23 


0.62 


0.70 


0.66 


e + e~ 


v^ = 4.16 GcV 


0.85 ±0.16 


0.63 


0.71 


0.66 


e + e~ 


= 5.20 GcV 


0.70 ±0.16 


0.59 


0.70 


0.28 


e+e~ 


Vs = 10.55 GeV 


0.56 ±0.08 


0.61 


0.74 


0.38 


e + e~ 


= 91.2 GeV 


0.60 ±0.04 


0.61 


0.75 


0.39 




E % = 200 350 GcV 


0.61 ±0.02 








7 A 




0.57 ±0.09 








i/iV 


E u ~ 25 GcV 


0.50 ±0.12 




0.69 


0.27 



Table 5: Summary of the available data on Fy and predictions from different models. The error 
for the theoretical predictions is not shown being relevant only for the third digit. 

For these reasons we assume that R\ measured in e + e~ (see Table |§|) can be used in neutrino induced 
charm-production, too. 

As an important by product of our study we have also extracted Fy from different processes and at 
several energies (see Table ||) . The simplest model to predict Fy is based on the spin-counting. Namely, 
vector mesons are spin-one states ( 3 5*i), while pseudo-scalar mesons are spin-zero states ( 1 S U ), therefore 
F v = 0.75. The discussion of more refined models (UCLA, JETSET, HERWIG and others) is beyond the 
purposes of this paper. For details we refer to ]17|] . 

From Table we can see that the measured Fy is independent of the processes and of the energy. This 
means that the probability for a c-quark to fragment into a D or a D* meson is universal and does not 
depend neither on the process nor on the energy. Notice that the UCLA model is the best in describing 
available e + e~ data. 

2.2 Measurement of the D s to D° and A c to D° ratios 

The ratio D s /D° has been measured in e + e~, ttN and 7 A?" experiments. A summary of the available data 
is given m Table |. 

The decay mode used to tag the event is 

D+ -» -> (K-K + )ir + 
whose BR, as reported by the Particle Data Group, are 

BR(D+ -> 07T+) = 0.036 ± 0.009 
BR(4> -» K~K + ) = 0.492 ± 0.007 

From Table || we can see that the D s /D° ratio is, within the errors, independent of the energy and of 
the process. 

Data on the ratio A c / D° are very poor. Indeed, it has been measured only in e + e~ experiments (see 
Table §). 
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(GcV) 


R s = D s /D" 


R c = A c /D { > 


4.14 


:] 


0.176 ±0.076 




10.55 




0.148 ±0.052 


0.148 ±0.052 


10.55 




0.184 ±0.040 


0.158 ±0.039 


M z o (ALEPH) 


0.206 ±0.039 


0.140 ± 0.017 


M z o (DELPHI) 


0.228 ±0.033 


0.158 ±0.039 


M z o (OPAL) 


0.153 ±0.047 


0.082 ±0.041 


WA92: 350 GeV/c ir- 
on Cu 


0.160 ±0.037 




WA92: 350 GeV/c ir- 
on W 


0.183 ±0.068 




NA14/2 (7JV) 


0.185 ±0.083 




Weighted average 


0.193 ±0.016 


0.135 ±0.015 



Table 6: Summary of the available R s and i? L 
experiments. The irN data have been taken from Refs. 
results from Ref. [Si. 



measurements extracted from e + e , z^N and 7N 
and references therein, while jN 



11 



The decay mode used to tag the event is 



A c — > pK 7T + 

whose BR, as reported by the Particle Data Group, is j|: 



BR(A C -» P K~ir + ) = 0.050 ± 0.013 

From Table | we can see that, although with a smaller statistical accuracy, both R s and R c are, within 
the errors, independent of the energy. 



2.2.1 Summary and discussion of all available data on R s and R c 

Although from Table || it seems that R s is constant over a wide range of energies (y/s ~ 4 4- 90 GeV) and 
independent of the process, some comments are in order. 

Given the quark composition of the Dj(cs) meson, it has to be created always together with at least 
one K meson. Therefore, being mj{ ~ 500 MeV, we expect that the threshold effect for D s production is 
more pronounced than for D mesons, when at least one 7T has to be produced (m n ~ 100 MeV). To account 
for the different threshold effect at low energies, in the following we do not use the R s values measured at 
the Z° peak. Furthermore, under the assumption that y/s in collider experiments can be replaced with W 
in fixed target experiments, W being the final state hadronic mass, and noting that neutrino-induced charm 
events at present experiments are characterized by values of W in the range 4 4 10 GeV, we can argue a R s 
value of R s = 0.171 ± 0.029 which corresponds to the weighted average of the first three results in Table ||. 

The available measurements on R c are very poor. Nevertheless, as it will be discussed in Section ^, we 
do not use the R c value to predict the charmed fractions in events induced by neutrinos. 
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3 Charm-production fractions in neutrino experiments 
3.1 The method 

From the previous Sections we can argue that once a charm-quark has been produced in deep-inelastic 
interactions (i.e. the energy of the process is higher than the threshold), it has a probability to produce a 
charmed hadron Ch which is, within the experimental errors, independent of the process and of the energy. 
Therefore, as far as the deep-inelastic scattering (DIS) is concerned, we can write the charm production rate 
as 



a c (DIS) _ a D o a D + a Ds cr Ac 



&CC &CC &CC &CC &CC 



°cc 



x(l + R 1 + R s + R c ) . (5) 



In the case of neutrino induced charm-production, Eq. (^|) is not correct. Indeed, in this case we also have 
to account for diffractive and quasi-elastic charm-production. Therefore, the inclusive charm-production rate 
in neutrino interactions can be written as 

gcQQ cr c (DIS) , 0" D a Ac 

= 1 \diff H \QE ■ (o) 

&CC ace a CC &CC 

If one accepts that Ri, R s and R c from e + e~ (with y/s = 4.1 4- 90 GeV) and other experiments can 
be used to described the fragmentation of charm-quarks produced in DIS neutrino interactions with average 
final state hadronic mass (W) ~ 10 GeV, then 



G ciy) _ <r D ° n p R p v °D S , °c 1 (7 s 

— X (1 + Ki + K s + K c ) H \diff H \qe {<) 

&CC &CC &CC &CC 

From Eq. (|7]) charm-production fractions in neutrino interactions can be written as 

f D o{E v ) = ^{E v )x 



&CC v ' 



f D+ (E u ) = R 1 x^(E v )x 



<Trr< V ' 



&CC 



f D+ (E u ) = (r s x ^(E u ) + ^ \ diff (E v )) 
V a cc &cc J 



<?cc y ' 



Notice that, given the poor knowledge on R c and on the quasi-elastic charm-production cross-section, 
we derive /a c by using the normalization constrain fpo + fjj+ + f D + + = 1. In order to estimate the 
charmed fractions and their energy dependence, we use the following inputs 

• the inclusive charm-production rate derived in Ref. p8[ ; 
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E v (GcV) 


Id° 


fo+ 


fnt 


/a+ 


5 


0.32 ±0.05 


0.12 ±0.02 


0.054 ±0.015 


0.50 ±0.18 


10 


0.46 ±0.07 


0.18 ±0.03 


0.078 ±0.022 


0.29 ±0.11 


15 


0.50 ±0.08 


0.20 ±0.03 


0.13 ±0.04 


0.18 ±0.07 


20 


0.52 ±0.09 


0.20 ±0.04 


0.13 ±0.04 


0.14 ±0.05 


25 


0.53 ±0.10 


0.21 ±0.04 


0.14 ±0.04 


0.12 ±0.05 


30 


0.54 ±0.10 


0.21 ±0.04 


0.14 ±0.04 


0.11 ±0.05 


35 


0.54 ±0.10 


0.21 ±0.04 


0.14 ±0.04 


0.11 ±0.04 


40 


0.54 ± 0.11 


0.21 ±0.04 


0.14 ±0.04 


0.11 ±0.04 


50 


0.55 ±0.12 


0.21 ±0.05 


0.14 ±0.04 


0.10 ±0.04 


60 


0.55 ±0.12 


0.22 ±0.05 


0.14 ±0.04 


0.09 ±0.04 


70 


0.56 ±0.13 


0.22 ±0.05 


0.14 ±0.04 


0.08 ±0.04 


80 


0.57 ±0.14 


0.22 ±0.06 


0.14 ±0.04 


0.06 ±0.03 


90 


0.58 ±0.15 


0.23 ±0.06 


0.15 ±0.04 


0.04 ±0.02 


100 


0.60 ±0.16 


0.23 ±0.06 


0.15 ±0.05 


0.02 ±0.01 


CHORUS 


0.524 ±0.036 


0.204 ±0.014 


0.128 ±0.012 


0.147 ±0.008 



Table 7: Prediction of charm-production fractions in neutrino induced events as a function of the 
neutrino energy. 



the energy dependence of the D° production rate reported in Ref. |L2| properly scaled to account for 



the effect discussed in Section 2.1.3; 



the energy dependence of the diffractive D s production given in Jl9[] , properly scaled in order to re- 
produce the average diffractive charm-production cross-section as measured by BEBC and NuTeV |2(| 
(g£ \diff= (0-31 ± 0.05) xlO- 2 ). 



3.2 Results on f h and and comparison with the data 

By using the method described in the previous section, we derived the charm-production fractions, as a 
function of the neutrino energy, as reported in Table and shown in Fig. [j]. Our results are in good agreement 
with the charm-production fractions extracted from the E531 data, see Fig. [|. It is worth noticing that + 
shows a dependence on the energy, higher values at low neutrino energies, consistent with the expectations. 
Indeed, quasi-elastic charm-production, which yields only A c , is expected to largely contribute to o c {y) for 
E v < 25 GeV g] . In Table ^ the expected charm- production fractions in the CHORUS experiment are also 
given. 

Having determined the /^'s, we can also estimate the semi-muonic branching ratio B^ of the charmed 
hadrons as a function of the neutrino energy. B^ is a very important quantity, being the input variable 
needed to extract from the dimuon data the element of the CKM matrix V c d- 

Recently, a direct measurement of B^ has been performed by the CHORUS Collaboration by using a 
statistics of about 1000 charm events reconstructed in the nuclear emulsions. Out of these, (88 ± 10 ± 8) 
dimuon events have been reconstructed, which correspond to [|13| 



B^ = (9.3 ± 1.3)% . 
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Figure 1: Predicted charm-production fractions as a function of the neutrino energy. The continous line 
shows the central value of the fh as measured by E531, while the dashed lines the upper and lower bounds 
at 68% CL. 
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Collaboration 


E v (GcV) 


B„ | Vcd | 2 




CDHS |2lTj 


20.0 


(0.41 ±0.07) x 10~ 2 


0.083 ±0.014 


CHARM II ]22] 


23.6 


(0.442 ±0.049) x 10~ 2 


0.090 ±0.010 


NOMAD m 


23.6 


(0.48 ±0.17) x 1Q- 2 


0.097 ±0.034 


CHORUS p[ 


27 




0.093 ±0.013 


CCFR (LO) |24 


140 


(0.509 ±0.036) x 10-* 


0.103 ±0.007 


CCFR (NLO) |25| 


140 


(0.534 ±0.060) x W- 2 


0.108 ±0.012 


PDG m 

(From unitarity at 90%) 


^ = 0.219-7-0.225 


{Yd) = 0.222 ±0.003 





Table 8: Charmed hadron semi-muonic branching ratios and Vcd measured by various experiments 
at different neutrino energies. The direct measurement performed by CHORUS is also shown. 

This measurement has to be compared with our prediction obtained by convoluting the charm-production 
fractions with the CHORUS neutrino flux 

= (8.8 ± 1.0)% . 

In Table ^ we derived, by using the V c( i value obtained by imposing the unitarity constraint to the CKM 
matrix and the measurements of B^ | V c d | 2 from various experiments, B^. Given the fact that the different 
experiments exploit different neutrino energy spectra, we can probe the sensitivity of B^ to the neutrino 
energy. As expected, the higher the neutrino beam energy the larger the value of B^. 

4 Conclusions 

We have presented a method to extract the charmed fractions in neutrino induced events. The method 
relies on the fact that, apart from processes peculiar to neutrinos such as quasi-elastic and diffractive charm- 
production, the charm-production and fragmentation mechanism is believed to be process-independent. We 
have verified this natural assumption going through a complete review of the available data from different 
experiments. Moreover, the D + over D° ratio is constant over the large energy range spanned by the collider 
and fixed target experiments reviewed in this paper. By using recent data from neutrino experiments, we 
have assessed the consistency of this ratio with the one predicted by other experiments. On the other side, 
the energy-independent behavior of the ratio itself is clear from the review of the experiments. 

Threshold effects for the D s production are seen and accounted for. By introducing the diffractive 
charm-production and using the unity constraint, we have predicted the charm fragmentation as a function 
of the neutrino energy in the range useful to present neutrino experiments. In particular, a prediction for 
the CHORUS experiment has been made. The determination of the fragmentation function also allows the 
prediction of the semi-muonic branching ratio of charmed hadrons. The prediction given is in good agreement 
with a recent measurement made by the CHORUS experiment. 
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